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fi 要 : Ak 
气象 数据 等 


于 遥感 与 CIS 技术 的 研究 进展 ,采集 涨 河流 域 高 程 数 据 \ 土 地 利用 


类 型 .土壤 类 型 特性 


在 SWAT 模型 中 模拟 与 面 源 污染 有 关 的 降水 .地 表 径 流 、 壤 中 流 、 下 渗 和 蒸 散 发 等 水 


MALE ,计算 侵蚀 量 .总 所 和 总 磷 等 面 源 污染 负荷 ,研究 涨 河流 域 面 源 污 染 特征 。 结 果 表 明 : 涨 河流 


域 多 年 平均 侵蚀 量 为 0.70 x107 t, 多 年 平均 总 


源 污染 防治 对 策 及 建议 。 
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Micaela 控制 水 平 的 提高 , 面 源 污染 
题 逐渐 显现 ,成 为 世界 范围 水 体 的 主要 污染 
a 面 源 污 染 是 降雨 (或 融雪 ) 产 流 或 灌溉 过 程 
中 ,冲刷 .侵蚀 等 使 地 表土 壤 累 积 的 污染 物 进 入 河 
流 .湖泊 等 受 人 随机 
性 和 不 确定 性 等 特征 -1。 全 球 30% ~ 50% 地 表 
水 在 不 同 程度 上 受到 面 源 污染 影响 ,美国 .德国 . 丹 
BE .芬兰 .英国 日 本 和 俄罗斯 等 国家 的 研究 也 表明 , 
面 源 污染 是 导致 水 环境 日 益 恶 化 的 主要 因素 之 
5 9。 我 国 太湖 SRR) 滇池 面 源 污染 问题 十 分 严 
EU ,长 江 流域 .黄河 流域 ` 松 辽 流 域 和 珠江 流域 
的 面 源 污染 状况 也 不 容 忽 视 呈 -5 。 由 现 期 看 ,我 国 
AREOSA E DEN 
污染 的 影响 "91 。 忽 视 面 源 污染 控制 ,不 利于 重点 污 
a 制度 的 顺利 实施 ,也 难以 达到 生 
态 环境 规划 目标 。 
王 营 等 "研究 了 涨 河流 域 土地 利用 现状 及 覆 
盖 植 被 变化 的 水 文 过 程 响 应 ; 朱 丽 等 "研究 了 近 
50 a 来 涨 河流 域 气候 变化 和 干旱 演变 过 程 ,发现 近 
41 a 涨 河流 域 大 多 数 年 份 无 旱 , 仅 个 别 年 份 出 现 干 
旱 , 且 程度 较 轻 ,下 游 地 区 干旱 相对 较 重 。 景 春 
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“ 氮 面 源 污染 负荷 为 244. 11 t, 涨 河流 域 多 年 平均 总 磷 
面 源 污染 负荷 为 27.91t, 涨 河流 域 1967 一 2017 年 多 年 平均 侵蚀 量 、 总 氮 面 源 污染 
负荷 总 体 呈 不 明显 的 波动 减少 趋势 ,同时 提出 了 改变 土地 利用 方式 .设置 植被 过 


、 总 磷 面 源 污染 
ee 坡 改 梯 等 面 


刚 ”发 现 涨 河流 域 水 环境 持续 向 好 ,干流 水 质 均 能 
达到 目标 值 要 求 ; 部 分 支流 受到 轻 度 污染 ,超标 物 为 
氨氮 。 目 前 对 于 涨 河流 域 面 源 污 染 与 土壤 侵蚀 量 、 
总 磷 .总 氮 关 系 的 模拟 却 少 有 研究 。 本 研究 拟 以 涨 
Bd 采用 ArcGIS 地 理 信息 系统 平 
A ,划分 流域 边界 ,解析 流域 内 土地 利用 类 型 .土壤 
nae 性 和 地 形 特征 ;通过 建立 流域 SWAT 模 
KJI ,模拟 与 面 源 污染 有 关 的 降水 地表 径 流 、 壤 中 
流 等 水 文 过 程 ,计算 侵蚀 量 .总 所 和 总 磷 等 面 源 污染 
负荷 ,研究 涨 河流 域 面 污染 特征 ;同时 从 景观 管理 
和 耕作 管理 等 方面 提出 涨 河流 域 面 源 污染 防治 措 
施 。 本 研究 有 利于 加 深 认识 流域 面 源 污染 机 理 , 扩 
展 遥 感 和 地 理 信 息 系统 技术 在 流域 面 源 污染 定量 
究 的 应 用 范围 ,为 提升 流域 环境 质量 和 生态 健康 水 
平 提供 有 效 的 方法 和 途径 。 


1 沙河 流域 概况 及 数据 来 源 


1.1 涨 河流 域 概况 
涨 河 王 位 于 中 国 甘 肃 省 南部 ,是 黄河 上 游 第 一 
大 文 流 , 发 源 于 青海 省 海南 州 蒙古 族 自治 县 境内 的 
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等 地 ,在 临夏 州 永 靖 县 刘 家 峡 水 库 大 坝 上 游 汇 人 黄 
河 。 干流 河道 长 673 km, 流域 面积 25 527 km’, WW 
门 村 站 多 年 平均 流量 172 m + s ,年 径流 量 53 x 
10° m (红旗 水 文 站 )。 自 产 水 资源 总 量 47.2 x 10° 
m ,为 甘肃 中 部 提供 了 丰富 的 水 资源 。 
1.2 数据 来 源 

涨 河流 域 流域 高 程 数据 DEM 源 自 ASTERC- 
DEM 数据 集 ,该 数据 集 由 NASA 对 地 观测 卫星 
Terra 热 发 射 和 反射 辐射 仪 观测 结果 进行 加 工 得 来 ， 
是 全 球 空间 分 辩 率 为 30 m 的 数字 高 程 数据 产品 。 
涨 河流 域 1: 10 万 土地 利用 源 自 寒 区 旱 区 科学 数据 
中 心 的 WestDC Version 1.0 1 km 分 辨 率 土 地 履 盖 数 
据 集 。 涨 河流 域 土壤 数据 来 自 基 于 世界 土壤 数据 库 
(HWSD ) 的 中 国土 壤 数 据 集 。 在 中 国 气象 数据 网 查 
询 涨 河流 域内 且 有 历史 数据 的 气象 站 ,并 从 地 面 气 
候 资 料 日 值 数 据 集 提取 流域 上 游 玛 曲 气 象 站 
1967—2017 年 流域 中 游 距 县 气象 站 1951—2017 
年 以 及 下 游 临 涨 气 象 站 1951 一 2017 年 的 地 面 气候 
资料 ,包括 降水 气温、 相对 湿度 风速 和 日 照 等 。 


2 研究 方法 


本 研究 采用 流域 尺度 模型 SWAT 模型 描述 降 
雨 径流 ,土壤 侵蚀 和 污染 物 迁 移 等 过 程 ” ,在 具有 
不 同 土壤 类 型 .土地 利用 和 管理 条 件 特 征 的 大 尺度 
复杂 流域 一 洗 河 流域 内 ,模拟 土地 管理 措施 长 期 对 
产 水 , 产 沙 及 农业 化 学 污染 物 负 谷 等 的 影响 。 
SWAT 模型 中 植物 生长 模型 可 用 于 评估 根 区 水 分 和 
营养 物 的 迁移 ,植物 的 蒸 散 发 和 生物 量 / 产 量 。 各 水 
文 响应 单元 中 的 侵蚀 量 和 产 沙 量 通过 修正 的 通用 土 
坏 流失 方程 即 MUSLE 来 计算 。 最 终 模拟 流域 多 种 
形式 的 氮 / 倒 的 运 移 和 转化 ,党 养 物 进入 干流 河道 或 
者 通过 地 表 径流 和 壤 中 流向 下 游 迁移 结果 。 


3 结果 分 析 


3.1 流域 边界 

涨 河流 域 采 用 多 幅 遥 感 影像 图 合并 ,生成 流域 
边界 。SWAT 模型 以 整个 涨 河流 域 为 对 象 ,以 刘 家 
峡 水 电站 为 流域 出 口 ,进行 流域 划分 ,最 终生 成 子 流 
域 139 个 。 生 成 的 泥 河 流域 图 与 中 国 地 图 彩色 版 仁 
加 后 显示 出 了 整个 泥 河 流域 范围 (图 1) 。 涨 河流 域 


图 1 涨 河流 域 边界 示意 图 
Fig.1 The diagram of the Taohe River Basin Boundary 


面积 25 512.84 km’ ,边界 长 度 1 638. 12 km ,最 低 高 
程 1748 m, 最 高 高 程 4 866 m ,平均 高 程 3 049m。 
3.2 土壤 类 型 

沙河 流域 薄 层 土 所 占 面积 最 大 为 11 559. 92 
km 73) , 占 流域 总 面积 的 45. 27% ,主要 分 布 在 流域 
的 中 上 游 。 其 次 为 黑 钙 土 ,面积 为 5 754.91 km , 占 
流域 总 面积 的 22. 54% , 主要 分 布 在 流域 中 游 。 始 
成 土 面积 为 4 573. 36 km , 占 流 域 总 面积 的 
17.91% ,主要 分 布 在 流域 下 游 。 高 活性 淋 溶 土 面积 
为 1667.36 km , 占 流域 总 面积 的 6.53% 。 黑 土 面 
积 为 702.92 km , 占 流域 总 面积 的 2.75% . BEB 
面积 为 471. 84 km’ , 占 流域 总 面积 的 1.85% 。 潜 育 
土 面积 为 410.51 km , 占 流域 总 面积 的 1.61% 。 灰 
色 森 林 土 面积 为 237. 32 km’, 占 流域 总 面积 的 
0.93% 。 钙 积 土 面积 为 94. 40 km , 占 流域 总 面积 
的 0.37% 。 冲 积 土 面积 为 4. 69 km , 占 流域 总 面 
积 的 0.18%。 有 机 土 面 积 为 18.45 km , 占 流域 总 
面积 的 0.07% 。 涨 河流 域 土壤 类 型 占 比 见 图 2。 

涨 河流 域 土壤 类 型 对 应 于 SWAT 模型 ,生成 土 
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图 2 涨 河流 域 各 土壤 类 型 占 比 
Fig.2 Proportion of soil types in Taohe River Basin 
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壤 属 性 ,水 河 流域 土壤 / 定义 见 表 1 , 涨 河流 


J 星 
Tab.1 Definition of soil attributes in the 域 土壤 属性 见 表 2a 
Taohe River Basin 3.3 植被 覆盖 及 土地 利用 类 型 
变量 名 定义 > ` > > yE 
ERA EX 旱地 划分 为 耕地 ,编号 为 1。 有 林地 、 江 木林、 
HYDGRP 土壤 水 文 单元 (A、B、C 或 D),A 为 低 径流 潜力 地 和 其 他 林地 划分 地 ,编号 为 2 Tr BEE 
土壤 ,B 为 中 等 径流 潜力 土壤 ,C 为 较 高 径流 潜 ads 其 他 林 me | 分 为 林 be ° eden 
力 土壤 ,D 为 高 径流 潜力 土壤 度 草地 ,中 覆盖 度 草地 和 低 履 盖 度 草地 划分 为 草地 ， 
SOL-ZMX 土壤 剖面 的 最 大 根 底 深 度 / mm 编号 为 3。 河 渠 湖泊 \ 水 库 坑 塘 \ 永久 性 冰川 雪 地 
S E Pee AHJ R 人 > ` ` ii J 
satz JEJERE /mm 和 滩涂 划分 为 水 体 ,编号 为 4。 城 镇 用 地 .农村 居民 
SOL-BD 土壤 的 湿 容 重 ,为 固体 颗粒 的 质量 与 土壤 总 体 、 > 、 、 > 、 yh 
hae ii 点 和 其 他 建设 用 地 划分 为 建设 用 地 ,编号 为 5。 沙 
chen 4 沼泽 地 RE Jh kE Bt j 其 他 划分 》 
SOL-AWC 土屋 的 有 效 含水 量 , 等 = 间 持 水 量 减 去 永久 地 .沼泽 地 、 裸 土地 \ 裸 岩石 砾 地 和 其 他 划分 为 未 利 
凋 世系 数 /mm + mm 用 地 ,编号 为 6( 图 3)。 
SOL- RAH . h~! Walk Ty i 光 > 尖 
SOL-K 饱和 渗透 系数 人 mm h 将 涨 河流 域 土地 利用 顶 格 数据 转换 为 面 数 据 ， 
SOL-CBN 有 机 碳 含量 / % sity 5 : 
H Mp H 统计 各 种 土地 利用 类 型 的 面积 占 流域 面积 的 比例 。 
SOL-CLAY 粘 粒 含量 , 粒 径 <0.002 mm / % 
涨 河流 域 草 : 只 只 最 为 15 257. 35 
SOL-SILT 粉 粒 含量 , 粒 径 在 0.002 ~0.05 mm / % vem wi ante Ta BA Pr itr — A r 
a y z H N y N 
SOL-SAND ”砂粒 含量, 粒 径 在 0.05 ~2. 0mm / % km , 占 流域 总 面积 的 59. 76% ,在 整个 流域 分 布 的 
` y ` ` 2 y 
SOL-ROCK AURA, REE > 2.0mm / % 比较 均匀 。 其 次 为 林地 ,面积 为 5 604.77 km , 占 流 
SOL-ALB 湿 土 的 反照 率 域 总 面积 的 21. 97% , 主要 分 布 在 涨 河流 域 的 中 上 
USLE-K USLE 方程 中 的 土壤 侵蚀 K 因子 游 。 耕 地 面积 3 633.30 km , 占 流域 总 面积 的 
R2 ， 涨 河流 域 土壤 属 ! 
Tab.2 Soil attributes in the Taohe River Basin 
类 编 土壤 HYD SOL-Z1 SOL-72 SOL-BDI1 SOL-BD2 SOL-AWC  SOL-K1 SOL-K2 SOL-CBN1 SOL-CBN2 
号 名 称 GRP / mm / mm /g-com™ ]/g cm~? / mm: mm™! /mm-h™! /mm-h7! /% / % 
1 黑 钙 土 B 300.00 1 000.00 1.40 1.40 0.14 16. 00 63.40 1.61 0.73 
2 E B 300.00 1 000.00 1.40 1.40 0.15 186. 00 165.00 0.46 0.46 
3 WREE B 300.00 1 000.00 1.39 1.37 0.14 9.20 6.70 1.07 0.51 
4 冲积 土 C 300.00 1 000.00 1.41 1.41 0.15 74.00 22.00 0.60 0.40 
5 潜 育 土 C 300.00 1 000.00 1.33 1.31 0.15 61.00 11.00 1.86 0.73 
6 灰色 森林 土 B 300.00 1 000.00 1.37 1.35 0.15 41.00 123.00 1.60 0. 40 
7 WEE B 300.00 1 000.00 1.37 1.35 0.15 37.00 74.00 1.31 0.53 
8 WEE C 300.00 1 000.00 1.51 1.51 0.06 33.50 52.00 1.74 1.74 
9 高 活性 淋 深 土 B 300.00 1 000.00 1.41 1.36 0.15 62. 50 109. 00 0.63 0.33 
10 黑土 B 300.00 1 000.00 1.38 1.39 0.15 17.20 41.00 1. 66 0.58 
11 有 机 土 D 300.00 1 000.00 1.32 1.39 0.15 42.00 6.30 39. 40 38. 46 
类 编 土壤 SOL- ”SOL-SILTI SOL-SAND1 SOL-ROCK1 SOL-CLAY2 SOL-SILT2 SOL-SAND2 SOL-ROCK2 USLE-K1 USLE-K2 
号 名 称  CLAYI/ 和 /% / % / % / % / % / % / % 一 a 
1 REE 21.60 39. 40 39.00 5.20 26. 60 39.00 34. 40 15.40 0.14 0.17 
2 钙 积 土 21.00 38.00 41.00 7.00 23.00 37.00 40. 00 7.00 0.17 0.17 
3 RE 23.11 39.22 37.67 7.44 26.00 37.67 36. 33 17.67 0.16 0.17 
4 ”冲积 土 18.00 48.00 34. 00 15.00 18.00 46.00 36. 00 10.00 0.18 0. 18 
5 Wa 33.33 35.67 31.00 3.00 36.00 31.33 32.67 2.33 0.13 0.15 
6 灰色 森林 土 21.00 54.00 25. 00 10. 00 28.00 49.00 23. 00 10.00 0.16 0.19 
7 Ht 23.50 43.50 33. 00 8.00 22.00 45.00 33. 00 5.00 0.15 0. 18 
8 WEE 13.0 38.75 48. 25 14. 00 13.00 38.75 48. 25 14.00 0.15 0.18 
9 ”高 活性 淋 深 土 、22. 50 30.50 47.00 4.50 29.50 28.50 42.00 6.50 0. 16 0. 16 
10 黑土 22. 50 41.00 36. 50 .50 22.00 40. 50 37.50 4.50 0. 14 0. 17 
11 有 机 土 32. 00 33.00 35. 00 28.00 25.00 23.00 52. 00 5.00 0.12 0.12 
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Fig.3 Vegetation cover and land use type map 


in the Weihe River Basin 


14.24% ,主要 分 布 在 涨 河流 域 的 中 下 游 。 未 利用 地 
为 737.42 km , 占 流域 总 面积 的 2.89% 。 建 筑 用 地 
203.55 km , 占 流域 总 面积 的 0. 80% 。 水 体 86. 35 
km’ , 占 流域 总 面积 的 0.34% 。 
3.4 水 文 响应 单元 

沙河 流域 坡度 分 为 <2° (9% ) 2° (9%) ~ 6° 
(11% ) .6"(11% ) ~ 1$"(27% ) .1$"(27% ) ~ 25° 
(47% )、>25"(47% ) 五 个 坡度 级 。 流 域 土地 利用 、 
土壤 和 坡度 的 数据 库 构 建 完 毕 以 后 ,将 每 一 个 子 流 
域 的 土地 利用 .土壤 和 坡度 进行 组 合 , 创 建 独一无二 
的 水 文 响 应 单元 。 
3.5 面 源 污 染 负 荷 
3.5.1 水 文 过程 计 算 结果 及 分 析 将 涨 河流 域 的 
土地 利用 土壤 .气象 等 数据 导入 SWAT 模型 中 ,从 
1967 年 1 月 1 日 一 2017 年 12 月 31 日 ,在 计算 时 间 
步 长 内 贡献 于 河流 的 地 表 答 流量 、 壤 中 流量 及 流域 

(1) 降水 量 

降水 量 是 在 计算 时 间 步 长 内 ,由 流域 内 雨量 站 
网 的 实测 降雨 量 通 过 泰 森 多 边 形 法 计算 得 到 的 流域 
平均 降水 量 。 涨 河流 域 1967—2017 年 多 年 平均 降 
水 量 为 582. 8 mm ,年 降水 量 最 大 值 出 现在 1967 年 ， 
为 813.8 mm, 年 降水 量 最 小 值 出 现在 1997 年 ,为 
415.5 mm。 涨 河流 域 1967—2017 年 的 降水 量变 化 
如 下 图 4 所 示 。 

对 涨 河 流域 年 降水 量 进 行 非 参 数 Man-Kendall 
趋势 检验 52 ,通过 计算 可 知 ,S = -85,Z = -0.68 < 


0 ,但 未 通过 信和 度 90% (+ 上 1.28) .95% (+ 上 1.64) 和 
99% ( £2.32) 显著 性 检验 ,表明 泪 河 流域 1967 一 
2017 年 降水 量 呈 减少 趋势 ,但 这 种 趋势 不 明显 。 

(2) 地 表 径 流 贡 献 量 

SWAT 模拟 结果 中 涨 河 流域 1967—2017 年 降 
水 对 于 地 表 径 流 贡献 量变 化 如 下 图 5 所 示 。 涨 河流 
域 多 年 平均 地 表 径 流 贡 献 量 为 48.7 mm ,年 地 表 径 
流 贡 献 量 最 大 值 出 现在 1967 年 ,为 97.5 mm, 年 降 
水 量 最 小 值 出 现在 1987 年 ,为 18.0 mm, 

对 涨 河流 域 年 地 表 径 流量 进行 非 参 数 Man- 
Kendall 趋势 检验 ,通过 计算 可 知 ,S =0,2 =0, 但 未 
通过 信和 度 90% ( +1.28)、95% ( £1. 64) Fil 99% 
(+ 上 2.32) 显 著 性 检验 ,表明 涨 河流 域 1967—2017 
年 地 表 径 流量 呈 不 明显 的 波动 减少 趋势 。 

(3) 壤 中 流量 

壤 中 流 是 在 时 间 步 长 内 根 区 和 包 气 带 中 贡献 于 
河流 的 横向 流动 水 量 。 通 过 SWAT 模拟 , 涨 河流 域 
多 年 平均 壤 中 流量 为 103.1 mm, 壤 中 流量 最 大 值 出 
现在 1967 年 ,为 159.7 mm ,二 中 流量 最 小 值 出 现在 
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Fig.5 Changes of surface runoff in Taohe River Basin 


from 1967 to 2017 
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1997 年 ,为 52.9 mm。 涨 河流 域 1967 一 2017 年 的 壤 
中 流量 变化 如 下 图 6 所 示 。 

对 旋 河 流域 年 壤 中 流量 进行 非 参 数 Man-Ken- 
dall 趋势 检验 ,通过 计算 可 知 ,S = - 185, Z = 
-1.49 <0,12Z1 >1.28, 通 过 信和 度 90% (+ 上 1.28) 显 
著 性 检验 ,表明 涨 河流 域 1967—2017 年 壤 中 流量 总 
体 呈 减少 趋势 。 

(4) 下 渗 量 

下 渗 量 是 在 时 间 步 长 内 下 渗 到 根 区 的 水 量 。 通 
过 SWAT 模拟, 涨 河流 域 多 年 平均 下 渗 量 为 109. 8 
mm ,下 渗 量 最 大 值 出 现在 1967 年 ,为 193.2 mm, F 
渗 量 最 小 值 出 现在 1997 年 ,为 33.4 mm。 涨 河流 域 
多 年 下 渗 量 变化 如 下 图 7 所 示 。 

对 涨 河 流域 下 渗 量 进行 非 参 数 Man-Kendall 趋 
势 检 验 , 通 过 计算 可 知 ,S = -217,Z = -1.75 <0, 
IZI >1.64, 通 过 信和 度 95% ( £1.64) 显著 性 检验 , 表 
明 涨 河流 域 1967—2017 年 下 渗 量 总 体 呈 波动 减少 
趋势 。 

(5) ZERE 
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图 6 1967—2017 年 涨 河流 域 壤 中 流量 变化 


Fig.6 Changes of flow in soil in Taohe River Basin 
from 1967 to 2017 
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图 7 1967—2017 年 涨 河流 域 下 渗 量 变化 
Fig.7 Changes of infiltration in Taohe River Basin 


from 1967 to 2017 
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图 8 1967—2017 年 涨 河流 域 疗 散发 量变 化 


Fig.8 Changes of evapotranspiration in Taohe River Basin 


from 1967 to 2017 


量 和 植物 蒸腾 量 。 将 SWAT 模拟 的 涨 河流 域 
1967—2017 年 蒸 散发 量 导 和 统计 分 析 软 件 SPSS 进 
行 数据 分 析 可 得 所 ] , 涨 河流 域 多 年 平均 蒸 散 发 量 为 
320.9 mm, 薰 散发 量 最 大 值 出 现在 2012 年 ,为 
365.4 mm; 葵 散 发 量 最 小 值 出 现在 1981 年 ,为 
291.9 mm。 涨 河流 域 多 年 莹 散发 量变 化 如 下 图 8 
所 示 。 

对 涨 河流 域 蒸 散发 量 进 行 非 参数 Man-Kendall 
趋势 检验 ,通过 计算 可 知 ,S =393,Z =3. 18 > 0， 
IZI >2.32 ,通过 信和 度 99% ( +2.32) 显 著 性 检验 , 表 
明 涨 河流 域 1967 一 2017 年 蒸 散发 量 总 体 呈 波动 上 
升 趋势 。 

3.5.2 面 源 污染 负荷 计算 结果 及 分 析 

(1) 侵蚀 量 

各 水 文 响应 单元 内 的 侵蚀 量 通 过 MUSLE( Mod- 
ified Universal Soil Loss Equation) 计算 。USLE 中 将 
降水 量 作为 侵蚀 能 力 的 指标 , 而 MUSLE 中 运用 径 
流量 模拟 侵蚀 量 。 将 涨 河流 域 1967—2017 年 侵蚀 
量 ” 模 拟 结果 导入 统计 分 析 软 件 SPSS 进行 数据 分 
析 可 得 , 涨 河流 域 多 年 平均 侵蚀 量 为 0.70 x 10 t, 
侵蚀 量 最 大 值 出 现在 1967 年 ,为 1.35 x 10° i; 侵蚀 
量 最 小 值 出 现在 1997 年 ,为 0.25 x 10° t。 涨 河流 域 
多 年 侵蚀 量变 化 如 图 9 所 示 。 

对 涨 河流 域 年 侵蚀 量 进 行 非 参 数 Man-Kendall 
趋势 检验 ,通过 计算 可 知 ,S= -43,Z = -0.34 <0, 
IZI <1. 28, 未 通过 信和 度 90% ( +1. 28)、95% 
(+1.64) 和 99% (+2.32) 显 车 性 检验 ,表明 涨 河流 
Ja 1967 一 2017 年 侵蚀 量 总 体 呈 不 明显 的 波动 减少 
趋势 。 

(2) BA 

SWAT 模型 计算 的 总 氮 面 源 污染 负荷 是 在 计算 
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图 9 1967—2017 年 泪 河 流域 侵蚀 量变 化 


Fig.9 Changes of soil erosion in Taohe River Basin 


from 1967 to 2017 


步 长 内 ,结合 土壤 类 型 ,模拟 从 各 水 文 响应 单元 通过 
地 表 径 流 进 入 到 河道 中 的 总 氮 负 和 荷 量 。 将 涨 河流 域 
1967—2017 年 总 氮 面 源 污染 负荷 导入 统计 分 析 软 
件 SPSS 进行 数据 分 析 可 得 , 涨 河流 域 多 年 平均 总 氮 
面 源 污染 负荷 为 244. 11 t, 总 氮 面 源 污染 负荷 最 大 
值 出 现在 2012 年 ,为 1 122.51 t; 总 拓 面 源 污染 负 和 荷 
最 小 值 出 现在 1997 年 ,为 25.51 t。 涨 河流 域 多 年 
总 氮 面 源 污染 负荷 的 变化 如 下 图 10 所 示 。 

对 涨 河流 域 年 总 氮 面 源 污 染 负荷 进行 非 参 数 
Man-Kendall 趋势 检验 ,通过 计算 可 知 ,$ = - 43 ,Z = 
-0.34 <0,121 <1.28 ,未 通过 信和 度 90% ( +1.28), 
95% ( +1.64) 和 99% ( + 上 2.32) 显 著 性 检验 ,表明 涨 
河流 域 1967—2017 年 总 氮 面 源 污 染 负荷 总 体 呈 不 
明显 的 波动 减少 趋势 。 

(3) 总 磷 

SWAT 模型 计算 的 总 磷 面 源 污染 负荷 是 在 计算 
步 长 内 ,结合 土壤 类 型 ,模拟 从 各 水 文 响应 单元 通过 
地 表 径 流 进 入 到 河道 中 的 总 磷 负 和 荷 量 。 将 涨 河流 域 
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图 10 1967—2017 EIKT MRE AUR 
Fig. 10 Changes of total nitrogen in Taohe River Basin 
from 1967 to 2017 


1967—2017 4F X WE EL WATS Ae Ha gay SA SET oP OT AK 
件 SPSS HEAT AGED AT HI FS, DEVAL Tat ah as AEP E 
面 源 污染 负荷 为 27.91 t, 总 磷 面 源 污染 负荷 最 大 值 
出 现在 2012 年 ,为 61.23 ;总 磷 面 源 污染 负 丛 最 小 
值 出 现在 1997 年 ,为 5.1 t。 涨 河流 域 多 年 总 磷 面 
源 污染 负 奏 的 变化 如 下 图 11 所 示 。 
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图 11 1967—2017 年 泪 河 流域 总 磷 变 化 
Changes of total phosphorus in Taohe River Basin 


from 1967 to 2017 


Fig. 11 


对 涨 河流 域 年 总 磷 面 源 污 染 负荷 进行 非 参数 
Man-Kendall 趋势 检验 ,通过 计算 可 知 ,S= -61,Z = 
-0.49 <0,121 <1.28 ,未 通过 信和 度 90% ( + 上 1.28)、 
95% (+ 上 1.64) 和 99% (+ 上 2.32) 显 著 性 检验 ,表明 涨 
河流 域 1967 一 2017 年 总 磷 面 源 污 染 负 荷 总 体 呈 不 
明显 的 波动 减少 趋势 。 


4 结论 


(1) 使 用 SWAT 模型 计算 了 涨 河 流域 1967 一 
2017 年 水 文 过 程 ,包括 降水 量 、 地 表 径 流 贡 献 量 、 壤 
HRE, FAE A E, (EH Man-Kendall 检验 
法 对 其 进行 了 趋势 分 析 与 突变 点 分 析 。 结 果 表 明 : 
涨 河流 域 1967 一 2017 年 多 年 平均 降水 量 为 582. 8 
mm ,年 降水 量 最 大 值 出 现在 1967 年 ,为 813.8 mm, 
年 降水 量 最 小 值 出 现在 1997 年 ,为 415.5 mm, Yk 
河流 域 多 年 平均 地 表 径 流 贡 献 量 为 48.7 mm ,年 地 
表 径 流 贡 献 量 最 大 值 出 现在 1967 年 ,为 97.5 mm, 
年 降水 量 最 小 值 出 现在 1987 年 ,为 18.0 mm。 通 过 
SWAT 模拟 , 涨 河流 域 多 年 平均 壤 中 流量 为 103.1 
mm , 培 中 流量 最 大 值 出 现在 1967 年 ,为 159.7 mm, 
培 中 流量 最 小 值 出 现在 1997 年 ,为 52.9 mm。 涨 河 
流域 多 年 平均 下 渗 量 为 109.8 mm ,下 渗 量 最 大 值 出 
现在 1967 年 ,为 193. 2 mm, 下 渗 量 最 小 值 出 现在 
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1997 年 ,为 33.4 mm。 涨 河流 域 多 年 平均 蒸 散 发 量 
为 320.9 mm , 蒸 散 发 量 最 大 值 出 现在 2012 年 ,为 
365.4 mm; 蒸 散发 量 最 小 值 出 现在 1981 年 ,为 
291.9mm。 涨 河流 域 多 年 平均 下 渗 量 为 109.8 mm, 
下 渗 量 最 大 值 出 现在 1967 年 ,为 193.2 mm, 下 渗 量 
最 小 值 出 现在 1997 年 ,为 33.4 mm。 涨 河流 域 多 年 
平均 荣 散 发 量 为 320.9 mm, ZPR a OB LE 
2012 年 ,为 365.4 mm; 茸 散发 量 最 小 值 出 现在 1981 
年 ,为 291.9 mm, 

(2) EH SWAT 模型 计算 了 涨 河流 域 1967 一 
2017 年 侵蚀 量 .总 氮 面 源 污染 负荷 .总 磷 面 源 污 染 
负荷 ,总体 呈 波动 减少 趋势 。 结 果 表明 涨 河流 域 多 
年 平均 侵蚀 量 为 0.70 x10° t ,侵蚀 量 最 大 值 出 现在 
1967 年 ,为 1.35 x105t; 侵 蚀 量 最 小 值 出 现在 1997 
年 ,为 0.25 x 10° ts。 涨 河流 域 多 年 平均 总 氮 面 源 污 
染 负荷 为 244. 11 it, 总 氮 面 源 污染 负荷 最 大 值 出 现 
在 2012 年 ,为 1 122. 51 t; MAUR TS Be it fi eM 
出 现在 1997 年 ,为 25.51 t。 涨 河流 域 多 年 平均 总 
磷 面 源 污 染 负 荷 为 27.91 t ,总 磷 面 源 污染 负荷 最 大 
值 出 现在 2012 年 ,为 61. 23 ti; 总 磷 面 源 污 染 负 荷 最 
小 值 出 现在 1997 年 ,为 5.1 to 


5 面 源 污染 防治 对 策 


为 进一步 降低 面 源 污 染 中 土壤 侵蚀 程度 及 总 
Be .总 须 负 和 荷 ,结合 涨 河 流域 特点 ,提出 以 下 对 策 : 

(1) 改变 土地 利用 方式 

将 易 造 成 水 土 流失 的 耕地 因地制宜 地 植树 造 
BK ,恢复 为 森林 植被 。 退 耕 还 林 之 后 地 表 径 流量 会 
相应 减少 ,流域 内 土壤 水 分 持 蕾 能力 的 提高 ,可 以 有 
效 地 减少 土壤 侵蚀 ,从 而 降低 泥 沙 的 输出 负荷 。 随 
着 耕地 面积 的 减少 ,流域 内 化 肥 使 用 量 将 会 大 幅度 
降低 , 氮 磷 营养 物 的 输出 也 会 大 幅度 的 前 减 。 

(2) 设置 植被 过 滤 带 

植被 过 滤 带 是 植物 密集 条 带 ,用 于 拦截 上 坡 污 
染 源 处 的 径流 ,并 进行 过 滤 。 植 被 过 滤 带 通过 降低 
坡 面 地 表 径流 的 流速 进而 导致 微粒 沉降 来 去 除 污 染 
物 。 同 时 ,植被 过 滤 带 也 可 以 增加 区 域 的 下 渗 ,减少 
地 表 径 流 和 非 颗粒 污染 物 。 植 被 过 滤 带 的 覆盖 植物 
常 选用 草 和 灌木 ,一般 选 址 于 河流 湖 库 岸 边 ,道路 、 
农田 的 周边 ,作为 坡 面 径流 的 进入 水 体 前 的 缓冲 。 

(3) 坡 改 梯 

梯田 是 田间 用 于 拦截 径流 和 预防 侵蚀 的 堤坝 。 
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坡 改 梯 是 在 坡地 上 分 段 沿 着 等 高 线 建造 阶梯 式 农 
H ,以 减少 片 蚀 和 沟 蚀 作用 ,增加 水 保 效 果 。 
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SWAT simulation and countermeasures of non-point source 


pollution based on remote sensing and GIS in Taohe River Basin 


LIU Quan-yang, QI Ming-liang, MA Xiao-zhou, GONG Juan-xiao 
(Gansu Academy of Environmental Sciences , Lanzhou 730000 , Gansu , China) 


Abstract: Based on the research progress of remote sensing and GIS technology, the elevation data, land use 
types, soil type characteristics , and meteorological data in the Taohe River basin in Qinghai and Gansu Province, 
China were collected. The precipitation , surface runoff, soil flow , infiltration ,and evapotranspiration parameters rela- 
ting to non-point source pollutions were simulated using SWAT model to calculate non-point source pollution loads 
of the erosion, total nitrogen and total phosphorus, and study the non-point source pollution characteristics in the 
Taohe River basin. The hydrological cycle simulated by SWAT model was based on the water balance equation to 
predict the runoff of each hydrological response unit. We then calculated the total runoff of the basin, which im- 
proves the prediction accuracy and more accurately describes the physical processes of the water balance. The 
effects of long-term land management measures on water production , sediment yield , and agricultural chemical pollu- 
tant loads were predicted in the large-scale complex watershed of Taohe River basin with different types of soil, land 
uses ,and management conditions. The plant growth model in the SWAT model can be used to assess water and nu- 
trient migration , plant evapotranspiration ,and biomass/yield in the root zone. The intensity of erosion and sediment 
yield in each hydrological response unit were calculated using the modified universal soil loss equation-MUSLE. Fi- 
nally , various forms of nitrogen/phosphorus migration and transformation were simulated in the basin, with nutrients 
entering the mainstream and migrating downstream via surface runoff and subsurface flows. Studies showed that the 
average annual erosion was 0.70 x 10° t and the average annual total nitrogen and phosphorus of non-source pollu- 
tion loads were 244.11 t and 27.91 t,respectively in the Taohe River basin. The average annual erosion, total nitro- 
gen,and total phosphorus of non-point source pollution loads showed a decreasing trend with indistinct fluctuations. 
At the same time, countermeasures and suggestions for non-point source pollution prevention and control, such as 
changing land-use patterns, setting up vegetation filter belts ,and changing slopes to terraces were proposed. 


Key words: Taohe River Basin;non-point source pollution; SWAT model; countermeasures 


